Mesons and tachyons with confinement and chiral restoration, and NA60 
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In this paper the spectrum of quark-antiquark systems, including light mesons and tachyons, is 
studied in the true vacuum and in the chiral invariant vacuum. The mass gap equation for the 
vacua and the Salpeter-RPA equation for the mesons are solved for a simple chiral invariant and 
confining quark model. At T=0 and in the true vacuum, the scalar and pseudoscalar, or the vector 
and axial vector are not degenerate, and in the chiral limit, the pseudoscalar groundstates are 
Goldstone bosons. At T=0 the chiral invariant vacuum is an unstable vacuum, decaying through 
an infinite number of scalar and pseudoscalar tachyons. Nevertheless the axialvector and vector 
remain mesons, with real masses. To illustrate the chiral restoration, an arbitrary path between 
the two vacua is also studied. Different families of light-light and heavy-light mesons, sensitive to 
chiral restoration, are also studied. At higher temperatures the potential must be suppressed, and 
the chiral symmetry can be restored without tachyons, but then all mesons have small real masses. 
Implications for heavy-ion collisions, in particular for the recent vector meson spectra measured by 
the NA60 collaboration, are discussed. 



I. INTRODUCTION 

Very recently, the precise di-muon measurement in 
heavy ion indium- indium collisions by NA60 Q, 0, S 0, 
collaboration provided an exceptional probe to observe 
vector mesons in excited vacua. The masses of vector 
mesons in excited vacua, have been extensively modelled, 
with different results, since Brown and Rho |E S IE 
proposed the scaling of the light-light vector mesons with 
the restoration of chiral symmetry. Notice that tachyons 
may also occur. When there are only mesons in the vac- 
uum, the vacuum is a minimum. It is then stable when 
the minimum is absolute, or metastable when the mini- 
mum is local because then the vacuum can decay through 
tunnelling. However when both mesons and tachyons oc- 
cur, the vacuum is a saddle point. The tachyons indicate 
the decay directions of the vacuum, and thus the vac- 
uum is unstable, it is a false vacuum. Here I compare 
the mesons and tachyons in the false chiral invariant vac- 
uum and in the true vacuum, in the framework of a chiral 
invariant and confining potential. 

Notice that, in the true vacuum of QCD, quarks are 
confined. On the other hand, in the excited vacuum, 
chiral restoration is expected. Therefore a framework 
with a confining and chiral invariant quark interaction 
is convenient to study mesons and tachyons in the two 
vacua. The present study, with confinement, upgrades 
our knowledge of vacua and of vacuum fluctuations in 
hadronic models. For instance vacua properties of the 
non-confining sigma model 0, 0, an d Nambu and 
Jona-Lasinio model 01 have been explored in detail , 
including suprising unstabilities led by the 't Hooft 
breaking determinant |l5j . In the simplest scenarios, the 
vacua manifold of these models has the well known Mex- 
ican hat shape, where the chiral invariant unstable vac- 
uum has a finite number of tachyons. The tachyons in 
the flavour SU(2) sigma model occur in the scalar a and 
in the pseudoscalar 7r + , ir°, ir~ channels. So in the sigma 
model there are four tachyons in the false chiral invariant 
vacuum, while in the true chiral symmetry breaking vac- 



uum there is one massive meson, the scalar a and three 
pseudoscalar mesons 7r + , 7r° , u~ . In the chiral limit the 
pseudoscalar mesons are goldstone bosons, in the border- 
line between mesons and tachyons. 

However, when quarks suffer a confining potential, 
the tachyon structure of the false chiral invariant vac- 
uum possibly differs from the one of the sigma model or 
the one of the Nambu and Jona-Lasinio model. In the 
true vacuum, the confining quark models have an infi- 
nite number of states in each channel, while the sigma 
model of the Nambu and Jona-Lasinio model only have 
a finite number of mesons. Other differences also oc- 
cur. Le Yaouanc et al. ^(| found that, even at high 
temperatures, the confining potential prevents a phase 
transition from the chiral symmetry breaking vacuum to 
the chiral invariant vacuum. Le Yaouanc, Oliver, Ono, 
Pene and Raynal, 0] also found that, with a harmonic 
confinement, there is an infinite tower of excited vacua, 
interpolating between the true chiral symmetry breaking 
vacuum to the highest chiral invariant vacuum. This re- 
sult was recently generalized to any confining potential 
by PB and Nefediev 0] . The existence of tachyons in the 
chiral invariant vacuum of a confining quark model was 
already signalled by Le Yaouanc, Oliver, Ono, Pene and 
Raynal, |l7|. Here these tachyons are studied in detail. 

Because the present problem is quite technical, and 
because it is not clear yet what is the best chiral invari- 
ant and confining quark model, for clarity I now use the 
framework of the simplest confining and chiral invariant 
quark model OS [II • 

Notice that a calibration problem exists in chiral com- 
putations. The full hadron spectrum remains to be cor- 
rectly reproduced. When the quarks were discovered, 
it was realized that the main difficulty of the quark 
model consisted in understanding the low pion mass. 
But Nambu and Jona-Lasinio [pj had already shown 
that the spontaneous dynamical breaking of global chi- 
ral symmetry provides a mechanism for the generation of 
the constituent fermion mass and for the almost vanish- 
ing mass of the pion. This mechanism was extended to 
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the quark model by le Yaouanc, Oliver, Ono, Pene and 
Raynal with the Salpeter equations in Dirac structure 
[l7j and by PB and Ribeiro with the equivalent Salpeter 
equations in a form 0] identical to the Random Phase 
Approximation (RPA) equations of Llanes-Estrada and 
Cotanch [2]]]. These chiral quark models also comply 
with the PCAC theorems, say the Gell-Mann Oakes and 
Renner relation [HIES, the Adler Zero plIHI^. the 
Goldberger-Treiman Relation [22I |25| , or the Weinberg 
Theorem |22, H3, H(| . Possibly a chiral quark model with 
the correct spin-tensor potentials will eventually repro- 
duce the full spectrum of hadrons [l^ . Nevertheless this 
is only a quantitative problem, qualitatively the simple 
model used here is sufficient to study several implications 
of chiral symmetry and confinement. 

Recently, the full mesonic spin-tensor potentials of the 
present simple model were determined for a quark and 
an antiquark with different isospin |27| . Here I exactly 
solve these boundstate equations of mesons and tachyons 
in different vacua. Importantly, the hamiltonian of this 
model can be approximately derived from QCD, 



TABLE I: Matrix elements of the spin-dependent potentials 



H = I d x 



^(x) (m o -ia-V) i>{x) + -g 2 J d 4 y 

^)yY^www^(»))%)y^(y) + ••• (i) 



up to the first cumulant order, of two gluons [281129. 30, 
|3l|, which can be evaluated in the modified coordinate 
gauge, 

g 2 (A^x)A b u (y)) ~ -^<W„o<M, [if 3 (x - y) 2 - U] (2) 

and this is a simple density-density harmonic effective 
confining interaction, mo is the current mass of the 
quark. The infrared constant U confines the quarks but 
the meson spectrum is completely insensitive to it. The 
important parameter is the potential strength Kq, the 
only physical scale in the interaction. In the true chiral 
symmetry breaking vacuum Kq ~ 0.3 ± 0.05 GeV fits 
reasonably the hadron spectra. However in the chiral 
invariant vacuum the potential strength Kq is supposed 
to be greatly suppressed. For simplicity, I will consider 
a vanishing light quark mo and all physical results will 
scale only with the potential strength Kq. 

I now address the meson and tachyon spectrum in dif- 
ferent vacua. In Section II the quark mass gap equation 
and the bound state quark-antiquark equation are re- 
viewed. In Section III the mass gap and boundstate equa- 
tions are solved numerically and the spectrum is studied 
in an arbitrary interpolation between the true and the 
chiral invariant vacuum. In Section IV the tachyons solu- 
tions of the boundstate equation are analytically studied. 
These first studies are performed at vanishing tempera- 
ture. However in heavy ion collisions finite temperatures 
are reached, sufficient for a QCD phase transition. The 
conclusion is presented in Section V, including the esti- 
mation of temperature effects on the spectra. 
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II. T = MASS GAP AND BOUNDSTATE 
EQUATIONS 

The relativistic invariant Dirac-Feynman propagators 
|17| . can be decomposed in the quark and antiquark 
Bethe-Goldstone propagators |2fJ, close to the formalism 
of non-relativistic quark models, 



& Dirac (k ,k) 



jk — m + it 

i 

k - E(k) + ie 
i 

~-ko-E{k) 



u«(k) 



w.(k) 



1 + S 



1 + S 



-fc-(T75 U s (0) 



1-5 



k ■ (T75 



Vs(0) 



-icr 2 75<( k ) i 



(3) 



where S — sin(y) 



C = cos(ip) = 



and ip is a chiral angle. In the non condensed vacuum, 
if is equal to arctan ^ , but tp is not determined from 
the onset when chiral symmetry breaking occurs. In 
the physical vacuum, the constituent quark mass m c (fc), 
or the chiral angle ip(k) = arctan !2sA*) ; j s a variational 
function which is determined by the mass gap equation. 
Examples of solutions, for different light current quark 
masses mo, are depicted in Fig. ^ For simplicity in the 
remaining iof this paper mo = will be assumed, nev- 
ertheless the effect of a finite current quark mass can be 
estimated with a small increase of the dynamically gen- 
erated constituent quark mass m c . 

Then there are three equivalent methods to find the 
true and stable vacuum, where constituent quarks ac- 
quire the constituent mass. One method consists in as- 
suming a quark-antiquark Pq condensed vacuum, and in 
minimizing the vacuum energy density. A second method 
consists in rotating the quark and antiquark fields with 
a Bogoliubov-Valatin canonical transformation to diago- 
nalize the terms in the hamiltonian with two quark or an- 
tiquark second quantized fields. A third method consists 
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FIG. 1: The constituent quark masses m c (k), solutions of the 
mass gap equation, for different current quark masses mo. 
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in solving the Schwinger-Dyson equations for the propa- 
gators. Any of these methods lead to the same mass gap 
equation and to the quark dispersion relation. Here I re- 
place the propagator of eq. JSJ in the Schwinger-Dyson 
equation, 



dw' d 3 k' 
2tt (2tt) 3 

u{k>) s ,^{k') s , v(k') s ,^(k') s , 



= u\(k) \ kk • a + mof3 

E 



w'-E(k') + it -w'-E(k') + ie 
E(k) = ul(k) | kk ■ a + m /3 1 "" " 
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E 



u(k') s ,u^(k') s 
w' - E(k') + : 



v{k') s ,v\k') s 



E(k') 



u s (k), (4) 



where, with the simple density-density harmonic interac- 
tion , the integral of the potential is a laplacian and 
the mass gap equation and the quark energy are finally, 



Aip{k) = 2kS(k) - 2m C(k) 
E(k) = kC(k)+m a S(k) 



2S(k)C(k) 
v '(kf C(kf 



(5) 



k 2 



U 
~2 



Numerically, this equation is a non-linear ordinary differ- 
ential equation, ft can be solved with the Runge-Kutta 
and shooting method. Examples of solutions for the cur- 
rent quark mass m c (k) = kt&rup, for different current 
quark masses mo, are depicted in Fig. ^ 

The Salpeter-RPA equations for a meson (a colour sin- 
glet quark-antiquark bound state) can be derived from 



TABLE II: The positive and negative energy spin- 
independent, spin-spin, spin-orbit and tensor potentials are 
shown, for the simple density- density harmonic model of eq. 
©. v?'(fc), C(k) and G{k) = 1 - S(k) are all functions of the 
constituent quark(antiquark) mass. 
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the Lippman-Schwinger equations for a quark and an an- 
tiquark, or replacing the propagator of eq. © in the 
Bethe-Salpeter equation. In either way, one gets (2j| 



+M(P) - E(k 1 ) - E(k 2 ) 

-M(P) - E(h) - E(k 2 ) 
rfik' 

— V(k-k') [u(k[)^(k',P)vHk' 2 ) 
v(k' 1 )(p- t (k\P)u\k l 2 ) 



(6) 



k+i 



k 2 = k — -j and P is the total 



</>+(k,P) = 

4>-\k,P) = 
x (k,P) = 

where k\ = 

momentum of the meson. Notice that, solving for %, one 
gets the Salpeter equations of Yaouanc et al. (17| . 

The Salpeter-RPA equations of PB et al. [n| and 
of Llanes-Estrada et al. [2l| are obtained deriving the 
equation for the positive energy wavefunction <p + and for 
the negative energy wavefunction <f>~ . The relativistic 
equal time equations have the double of coupled equa- 
tions than the Schrodinger equation, although in many 
cases the negative energy components can be quite small. 
This results in four potentials V a @ respectively coupling 
v a = r<fi a to j/ 3 . The Pauli a matrices in the spinors 
of eq. 10 produce the spin-dependent potentials of 
Table CH 

Notice that both the pseudoscalar and scalar equations 
have a system with two equations. This is the minimal 
number of relativistic equal time equations. However the 
spin-dependent interactions couple an extra pair of equa- 
tions both in the vector and axialvector channels. While 
the coupling of the s-wave and the d-wave are standard in 
vectors, the coupling of the spin-singlet and spin-triplet 
in axialvectors only occurs if the quark and antiquark 
masses arc different, say in heavy-light systems. I now 
combine the algebraic matrix elements of Table^with the 
spin-dependent potentials of Table [H] to derive the full 
Salpeter-RPA radial boundstate equations (where the in- 
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frared U is dropped from now on). I get the J 1 



1 Sq pseudoscalar (P) equations, 



dk 2 



the J p = + , 3 Pq scalar (S) equations, 
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the J p = 1 , coupled 3 Si and 3 Di vector (V and V*) equations , 
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the J p = 1 + , coupled x Pi and 3 Pi axialvector (A and A*) equations 
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III. NUMERICAL SOLUTION OF THE MASS 
GAP AND BOUNDSTATE EQUATIONS AT T — 



In the light-light limit of m q = m q —* and tp — > 0, 
it is clear that eq. (J7J) and eq. (JSJ) become identical. 
They also possess takyonic solutions [l^. I n the- same 
limit, eq. © can be block diagonalized 0], an d each 
block, with mixed s-wave and d-wave, is identical one of 
the two independent blocks of eq. (|10|) . This checks that 
the chiral partners P-S and V, V*-A, A* are degenerate 
in the false chiral symmetric vacuum. 

Another interesting case is the heavy-light case where, 
say, the antiquark has a mass m q ~ mo q >> Kq, there 
are no Tachyons, and the negative energy components 



nearly vanish, like in non-relativistic quark models. In 
the infinite m q limit, Sq — > 1, and the antiquark spin is 
irrelevant, see Tabled complying with the Isgur-Wise 
heavy-quark symmetry. 

Notice that this model, like any chiral model, has the 
same number of meson states in the spectrum as the nor- 
mal quark model. The mass splittings can de related, as 
usual, to spin-tensor potentials. 

For the numerical solution, I change the sign of the 
second and fourth lines in eqs Q to (fTUf) and then I get 
a simple eigenvalue equation. The results are shown in 
Figs. El El LU an d El In Fig. El the pseudoscalar and 
scalar light quark and light antiquark meson masses are 
interpolated from the true spontaneously chiral symme- 
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FIG. 2: Light-light meson masses, in (a), pseudoscalar, in 
(b), scalar, when the light quark mass interpolates from the 
zero mass of the chiral invariant false vacuum to the solution 
ra c of the mass gap equation in the true vacuum. The dark 
curves correspond to mesonic real masses and the light curves 
correspond to tachyonic imaginary masses. 




0.2 0.4 0.6 0.8 1 1.2 1.4 

m/m c 




0.2 0.4 0.6 0.8 1 1.2 1.4 

m/m c 



FIG. 3: Light-light meson masses, in (a), vector and in (b), 
axial, when the light quark mass interpolates from the zero 
mass of the chiral invariant false vacuum to the solution m c of 
the mass gap equation in the true vacuum. The dark curves 
correspond to mesonic real masses and the light curves corre- 
spond to tachyonic imaginary masses. 
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try breaking vacuum to the false chiral restored vacuum. 
In Fig. |3 the vector and axial light quark and light anti- 
quark meson masses are interpolated from the true spon- 
taneously chiral symmetry breaking vacuum to the false 
chiral restored vacuum. In Fig. 01 the pseudoscalar and 
scalar light quark and heavy antiquark meson masses are 
interpolated from the true spontaneously chiral symme- 
try breaking vacuum to the false chiral restored vacuum. 
In Fig. [5J the vector and axial light quark and heavy 
antiquark meson masses are interpolated from the true 
spontaneously chiral symmetry breaking vacuum to the 
false chiral restored vacuum. 

A remarkable result of the numerical finite differ- 
ence solutions is that all studied pseudoscalar and 
scalar mesons, including all radial excitations, become 
tachyons, with arbitrarily large imaginary masses. This 
will be confirmed in the next Section IV. 



On the other hand, all the other mesons suffer small 
mass over potential strength M/Ko corrections from one 
vacuum to the other. Notice however that the potential 
strength K$ is expected to change significantly when the 
vacuum is changed. This will be discussed in detail in 
Section V. 

A remarkable feature of the vector meson groundstate 
in Fig. may be relevant for the p and u> mesons. Al- 
though the M/Kq corrections from one vacuum to the 
other are small, at small but non- vanishing quark mass 
the groundstate vector meson is a tachyon. This occurs 
just before the vector and axial vector are degenerate. 
Because the actual light current quark masses are small 
but non- vanishing, this will be addressed in Section V 
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FIG. 4: Heavy-light meson masses minus the infinitely heavy 
antiquark mass, in (a), pseudo-scalar, in (b), scalar, when the 
light quark mass interpolates from the zero mass of the chiral 
invariant false vacuum to the solution m c of the mass gap 
equation in the true vacuum. In this case there are no tachy- 
onic imaginary masses. 
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FIG. 5: Heavy- light meson masses minus the infinitely heavy 
antiquark mass, in (a), vector and in (b), axial, when the 
light quark mass interpolates from the zero mass of the chiral 
invariant false vacuum to the solution m c of the mass gap 
equation in the true vacuum. In this case there are no tachy- 
onic imaginary masses. 

M 

10 r 



o 



a 



0.2 0.4 0.6 

- n 



1 1.2 1.4 

m/m c 



M — niq 



0.2 0.4 0.6 



1 1.2 1.4 

m/m c 



IV. ANALYTICAL STUDY OF THE TACHYONS 

I now study in detail the properties of the eigenval- 
ues of the Salpeter or RPA equations. The boundstate 
equation can be decoupled, 



Vv- 
Vv- 



Mv+ 
-Mv- 



{H + V){v+ +v~) = M{v 
(H — V)(v + - v ) = M(V 

{H-V){H + V){v+ + v-) 
{H + V){H -V){v+ - v ) 



= Af 2 (^ 
= M 2 (^ 



(11) 
v-) ■ 



Thus we get a pair of eigenvalue equations, where H and 
V are hermitean, but (H-V)(H+V) and (H+V)(H-V) 



are not hermitean. Nevertheless M 2 can be prooved to 
be real, if one considers two different eigenvalues Mi 2 

(h - v){H + 10(^1+ + vr)= m 1 2 (i/ 1 + + vr) 

{v 2 + - v 2 -)HH - V){H + V)= {v 2 + - v 2 -)^M 2 2 * 



(Mi 2 - M 2 2 *){is 2 + - v 2 -)\ Vl + + 1/1-)= 
(M 2 2 - Af 1 2 *)(i/ 1 + - vi-)\u 2 + + V2 ) = 



(12) 



Notice that the orthonormalization condition 0, IT^. l20j 
of the Salpeter- RPA equation is, 



1 

-1 



= 5.. 



1,3 



(13) 



Thus, either the two eigenvectors are orthogonal or the 
squared eigenvalue M 2 is real. This shows that the so- 
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lutions of the boundstate equation can only have real or 
purely imaginary masses. While the real masses corre- 
spond to mesons, the imaginary masses correspond to 
tachyons. 

I now study in detail the solutions in the chiral in- 
variant vacuum and in the chiral limit, where both the 
current mass m and the constituent mass m c vanish. In 
general the boundstate equations decouple in two differ- 
ent equations, one for J > with, 



H = - 
V= ± 



dk- 



+ 2k 



and another for J > 1 with, 



H = — 

V- -2- 

v fe 2 



2k 



k 2 



k! 2 



(14) 



(15) 



,2 -i 

Notice that the different potentials — 4u? ■> 2fc , A- are 
bound from below and positive definite in the sense that 
all their eigenvalues are positive. However — -X is un- 
bound from below. Thus, in eq. O all terms H + V or 
H — V are positive definite and bound from below, ex- 
cept for the H — V of the J = pseudoscalar and scalar 
tachyons in eq. I|14fl . 

Notice that Fig. |21 suggests that all pseudoscalars 
and scalars become tachyons in the chiral invariant vac- 
uum. To confirm this suggestion of an infinite number 
of tachyons it is convenient to regularize the scalar and 
pseudoscalar equations, because the wave- functions are 
concentrated at extremely small distances. A very small 
quark mass m is assumed constant for simplicity, and the 
momentum and mass are rescaled, 



k/m 
Mm 2 



k' , 
M' 



(16) 



Notice that any finite solution M' in fact corresponds to 
an infinite mass M = M'/m 2 , and that a wave-function 
with a finite k' corresponds to a wave-function with in- 



finitesimal momentum k 



one gets for the pseudoscalar, 



Then, starting from eq. 



H + V = 
H-V = 



JL 

' dK 
JL 
' dk r 



k 2 + l 



jwtw 



(17) 




respectively positive definite and with negative eigenval- 
ues, and from eq. one gets for the scalar, 



(18) 



respectively positive definite and with negative eigenval- 
ues. An irrelevant term m 3 2k' / \J k' 2 + 1 is also present 
in the rescaled equations. 

The Bohr-Sommerfeld quantization condition can be 
used to count the number of negative eigenvalues of the 



TABLE III: Masses of the first angular and radial excitations 
of the different light-light tachyons and mesons in the chiral 
invariant false vacuum and in the chiral limit. Each column 
includes both positive and negative parity degenerate states, 
except for the pseudoscalar and scalar tachyonic states. No- 
tice that the tachyon masses are infinite and that they are 
regularized by an arbitrarily small quark mass m. The meson 
masses are separated in two different families with the same 
J because two different equations decouple for each J. 
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H — V pseudoscalar operator and of the H + V scalar op- 
erator. The leading term at high momentum, assuming 
the highest possible negative mass M' ~ 0, is, 



1 



1 



dk! 



(19) 



This shows that the number of tachyons in the pseu- 
doscalar and scalar channels are both infinite. 

This is confirmed by the numerical solution of the reg- 
ularized Salpeter equation. In Table IIIII we show the 
masses of the different light-light tachyons and mesons in 
the chiral invariant false vacuum and in the chiral limit. 

Notice that there is no pseudoscalar-scalar degeneracy 
in the rescaled equations since the equations are differ- 
ent, M's 7^ M 1 ps- Nevertheless both pseudoscalar and 
scalar tachyons have infinite imaginary masses and we 
get Ms = Mp — ooi. 



CONCLUSION, INCLUDING 
TEMPERATURE EFFECTS 



Assuming a confining potential, the mass M spectrum 
of mesons is studied in the true chiral symmetry break- 
ing vacuum and in the unstable vacuum where chiral 
symmetry restoration occurs. The only parameter is the 
strength Kq of the potential. Chiral models have the 
same number of meson states in the spectrum as the nor- 
mal quark model. The mass splittings can de related, as 
usual, to spin-tensor potentials. In the limit of vanishing 
constituent quark masses, all spin-dependent potentials 
are quite simple, proportional to Kq 3 jk 2 . 

In the chiral limit the mesons suffer small M/Kq 
changes from one vacuum to the other, except for the 
J = pseudoscalars and scalars. All the J = mesons, 
including all possible radial excitations, are transformed 
in tachyons with infinite imaginary masses, when the true 
vacuum is replaced by the chiral invariant vacuum. An 
detailed analytical proof and a precise numerical study 
of the tachyons are also presented here. 
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However, before moving to the conclusions, these beau- 
tiful mathematical results should be matched with our 
knowledge the deconfined phase of QCD. 

My first comment concerns the calibration problem 
of any chiral symmetric model. The Sigma Model, the 
Nambu and Jona-Lasinio model and Chiral Lagrangian 
estimations are not confining and thus are not expected 
to address correctly hadrons with spin, angular or ra- 
dial excitations. The present model is adequate to study 
the angular or radial excitations of hadrons, and in this 
sense it already upgrades previous estimations of the me- 
son spectra in the chiral restored vacuum. Nevertheless 
the present density-density interaction suffers from un- 
calibrated spin-tensor potentials. But I submit that the 
under development chiral invariant quark models with a 
confining funnel interaction |2lL l33l ] a vector interaction 
[2fl or long range scalar interactions [2E can 
be correctly calibrated. Nevertheless, for a qualitative 
study, the present density-density harmonic confining in- 
teraction should be sufficient, since PB and Nefediev [Isj 
have shown that this interaction has similar mass gap 
solutions to the other possible confining potentials in 
Coulomb gauge QCD. 

My second comment concerns the parameters of the 
present model. The potential strength Kq, the dominant 
scale of the present study, is expected to change from the 
ordinary QCD vacuum to the deconfined phase of QCD. 
This is quite important because the meson masses scale 
with Kq. 

Mys first conclusion concerns corrections due to the 
current quark mass. The light current quark mass is 
small but not vanishing. The u or d quarks correspond 
to and increase the chiral limit quark constituent mass by 
1% to 2% of m c , while the s quark amounts to increase 
the chiral limit quark constituent mass by up to 50%. 
For instance in the true vacuum the s constituent quark 
mass is of the order of 1.5 m c , while in the chiral restored 
vacuum the s constituent quark mass is of the order of 0.5 
m c . These simple factors are sufficient to estimate from 
the Figs. El 01 and [SJ the masses of the vectors p, lo or 
4>, or of the pseudoscalar and vector D and D s , relevant 
for the new di-muon measurements of NA60. In the light- 
light systems, with aaorrf quark and a u or d antiquark, 
the number and the imaginary mass of pseudoscalar and 
scalar tachyons are not infinite, nevertheless they are very 
large. 

Interestingly, in Fig. |3]the vector meson has real mass 
for zero quark masses, but for a small mass the vector 
meson is a tachyon. Thus it is possible that the p meson, 
or the lo meson, simply disappear in the chiral restored 
vacuum. Because the quark mass interval, where the 
vector meson is a tachyon, is quite small, it is plausible 
that the p meson and the lo meson may have a different 
tachyonic behaviour, although the present study cannot 
explore the differences between the p and the lo. Notice 
that the NA60 collaboration saw differences between the 
production rate of the p and the lo Q, 0, 0, 0, H , but this 
may also be due to p interactions with 7r at the periphery 



of the deconfined QCD bubble [37|. 

My second conclusion is that the chiral invariant vac- 
uum is too unstable to be reached, unless confinement 
is lost. This is clearly signalled by the infinite, (or very 
large) number of infinite (or very large) imaginary mass 
of tachyons in the pseudoscalar and scalar channels. This 
extreme unstability confirms a result of Le Yaouanc et al. 
[l(| , who studied the deconfinement transition, using the 
present confining potential, and concluded that the tran- 
sition does not occur for any finite temperature. There- 
fore a change in the potential must happen before the 
chiral restoration transition occurs. This also confirms 
the lattice QCD simulations initiated by Kogut, Wyld, 
Karsch and Sinclair |38|, |39|, |40j, and the Schwinger- 
Dyson calculations initiated by Bender, Blaschke, Kali- 
novsky and Roberts, 0,^3 who also found a restoration 
of chiral symmetry coincident with the loss of confine- 
ment at temperatures of the order of 150 MeV. 

The third conclusion of this paper is that all the meson 
masses are much smaller in the high temperature chiral 
invariant vacuum, than they are in the low temperature 
symmetry breaking vacuum. Notice for instance that the 
apparently constant vector and axialvector masses of Fig. 
|3]are proportional to the potential strength Kq, thus they 
decrease when the potential strength decreases. This is 
an educated conclusion, based on Lattice QCD simula- 
tions of the dependence of the confining potential with 
temperature and also with dynamical fermions. when 
confinement is lost, |40| at temperatures of the order 
of 150 MeV, the strength of the potential is also de- 
creased. These two effects are necessary for chiral sym- 
metry restoration. 

Assuming these two changes, both in shape and 
strength of the potential, the spectra computed in this 
paper can be reinterpreted. Assuming that confinement 
disappears, the infinite number of infinite imaginary mass 
tachyons go away. Moreover, a smaller strength of the 
potential is also necessary to remove any tachyon in the 
chiral symmetric vacuum. Then the chiral symmetric 
vacuum is the only and true vacuum. Notice that, for 
light quarks, the largely dominant scale, including the 
scale ruling the constituent quark mass, is the strength 
of the potential. All the spectra are proportional to the 
strength of the potential, see Figs|3 El El and [S] Then, 
with a much weaker potential, the masses and widths of 
any possible mesons are much smaller (except for the con- 
tribution of the heavy quark mass, say the c mass in D 
or D s mesons) than the masses of ordinary mesons listed 
by the Particle Data Group • Thus the vector mesons 
identified by the NA60 collaboration, with masses close 
to the ordinary masses, are not expected to be probed 
inside the deconfined phase of QCD, where all mesons, if 
any, are much lighter. 
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